Dynamical Coulomb blockade and spin-entangled electrons 
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We consider the production of mobile and nonlocal pairwise spin-entangled electrons from tunnel- 
ing of a BCS-superconductor (SC) to two normal Fermi liquid leads. The necessary mechanism to 
separate the two electrons coming from the same Cooper pair (spin-singlet) is achieved by coupling 
the SC to leads with a finite resistance. The resulting dynamical Coulomb blockade effect, which we 
describe phenomenologically in terms of an electromagnetic environment, is shown to be enhanced 
for tunneling of two spin-entangled electrons into the same lead compared to the process where the 
pair splits and each electron tunnels into a different lead. On the other hand in the pair-split pro- 
cess, the spatial correlation of a Cooper pair leads to a current suppression as a function of distance 
between the two tunnel junctions which is weaker for effectively lower dimensional SCs. 

PACS numbers: 74.45.+C, 74.50.+r, 73.23.Hk, 03.65.Ud 



Introduction.- The controled creation of nonlocal en- 
tanglement is crucial in quantum communication as well 
as in quantum computation tasks p). Quantum entan- 
glement is further interesting in its own right since it 
leads to a violation of Bell's inequality Q. Several solid 
state entanglers, a device that creates mobile and nonlo- 
cal pairwise entangled electrons, were proposed recently 
EJ [j S A particularly interesting quantity is 
the spin of the electron which was shown to be a promis- 
ing realization of a quantum bit [9( . A natural source of 
spin entanglement is provided by Cooper pairs in an s- 
wave superconductor (SC) , since the Cooper pairs are in 
a spin singlet state. Weakly coupling the SC to a normal 
region allows for (pair)-tunneling of Cooper pairs from 
the SC to normal leads and single particle tunneling is 
suppressed at low energies below the SC gap. Subse- 
quently, Coulomb interaction between the two electrons 
of a pair can be used to separate them spatially leading 
to nonlocality 10] . To mediate the necessary interaction 
entangler setups containing quantum dots or which ex- 
hibit Luttinger liquid correlations 0, Q (e.g. nanotubes 
in the metallic regime) were proposed recently. 
In this letter we show that a considerably simpler ex- 
perimental realization can be used to generate the nec- 
essary Coulomb interaction between the electrons of a 
pair. Indeed if the normal leads are resistive a dynamical 
Coulomb blockade (CB) effect is generated with the con- 
sequence that in a pair tunneling process into the same 
lead the second electron still experiences the Coulomb re- 
pulsion of the first one, which has not yet diffused away. 
Natural existing candidates with long spin decoherence 
lengths (~ 100 /itm ll[) for such a setup are e.g. semicon- 
ductor systems tunnel-coupled to a SC, as experimentally 
implemented in InAs [l^, InGaAs or GaAs/AlGaAs 
[lij . Recently, 2DEGs with a resistance per square up 
to almost the quantum resistance Rq = h/e 2 ~ 25.8 kfl 
could be achieved by depleting the 2DEG with a volt- 
age applied between a back gate and the 2DEG [l5| . 
In metallic normal NrCr leads of width ~ 100 nm and 
length ~ 10 fim, resistances of R = 22 — 24 kSl have been 



produced at low temperatures. Even larger resistances 
R = 200 — 250 kf2 have been measured in Cr leads ^(| . 
We use a phenomenological approach to describe charge 
dynamics in the electromagnetic circuit which is de- 
scribed in terms of normal-lead impedances and junc- 
tion capacitances, see Fig. 1. The subgap transport of 
a single SN- junction under the influence of an electro- 
magnetic environment has been studied in detail 1171 

EH 

In order to create nonlocal entangled states in the leads 
we have to go beyond previous work to investigate the 
physics of two tunnel junctions in parallel with two dis- 
tinct transport channels for singlets. A Cooper pair can 
tunnel as a whole into one lead, or the pair can split 
and the two electrons tunnel to separate leads, leading 
to a nonlocal entangled spin-pair in the leads [l^. ^ n 
the case where the pair splits we find that the CB effect 
provided by the electromagnetic environment is uncorre- 
cted for the two electron charges. In contrast, if the two 
electrons tunnel into the same lead we find a dynamical 
CB consistent with a charge q = 2e, where e ist the ele- 
mentary charge. Thus the CB effect is twice as large for 
the unsplit process which enhances the probability for a 
nonlocal (pair-split) process. On the other hand we show 
that the spatial correlations of a Cooper pair results in 
a suppression factor for tunneling via different junctions 
which is weaker for lower dimensional SCs. 
Setup and formalism.- The setup is sketched in Fig. 1. 
The SC is held at the (electro-)chemical potential us by a 
voltage source V. The two electrons of a Cooper pair can 
tunnel via two junctions placed at points r± and on 
the SC to two separate normal leads 1 and 2 which have 
resistances R\ and i?2, resp. They are kept at the same 
chemical potential /!/ so that a bias voltage ji = jis — Mi 1S 
applied between SC and leads |2JJ. The system Hamilto- 
nian decomposes into three parts H = H e + H env + Ht- 
Here H e = Hs + 2 n =i 2 ^in describes the electronic 
parts of the isolated subsystems consisting of the SC and 
Fermi liquid leads n = 1, 2, with H hl = J2 pa £ P cl pcr c npcr , 
where a = (1,1). The s-wave SC is described by the 
BCS-Hamiltonian H s - fJ-s^s = Ska E k italka with 
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(V)ev=^ 



We remark that in our setup the SC is held at the con- 
stant chemical potential /is by the voltage source, see Fig. 
1. Therefore the charge relaxation of a non-equilibrium 
charge on one of the capacitors described by QJ does not 
influence the charge dynamics of the other junction and, 
as a consequence, H env = Yl n =i 2^ e ™, n Electron 
tunneling through junctions 1, 2 located at points ri, Y2 
of the SC nearest to the leads 1,2 is described by the 
tunneling Hamiltonian Ht = 53n=i 2 ^Tn + h.c. where 



H 



Tn 



(2) 



FIG. 1: Entangler setup: A superconductor (SC) with chem- 
ical potential fis is tunnel-coupled (amplitude to) via two 
points ri and r-2 of the SC to two Fermi liquid leads 1,2 with 
resistance Ri t 2- The two leads are held at the same chemical 
potential fii such that a bias voltage /1 = /is — ^; is applied 
between the SC and the two leads. The tunnel-junctions 1,2 
have capacitances Clo- 



the quasiparticle spectrum Ek — {£% + A 2 ) 1 / 2 where 
6c = k 2 /2m — pis- The electron creation (cj^) and an- 
nihilation (cko-) operators are related to the quasiparti- 
cle operators by the Bogoliubov transformation cw-j- / n = 



where Uk and Vk are the usual 



Wfc7k( T /i) ± v kll HL/r) , 
BCS coherence factors. To describe resistance and dis- 
sipation in the normal leads we use a phenomenological 
approach [22|. where the electromagnetic fluctuations in 
the circuit (being bosonic excitations) due to electron- 
electron interaction and the lead resistances are modeled 
by a bath of harmonic oscillators which is linearly cou- 
pled to the charge fluctuation Q n of the junction capac- 
itor n (induced by the tunneling electron). This physics 
is described by [23, [ZQ 
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The phase <p n of junction n is the conjugate variable to 
the charge satisfying [4>„,Q m ] — ie8 n>m . As a conse- 
quence e~ i ^ n reduces Q n by one elemantary charge e. 



Here to is the bare electron tunneling amplitude which 
we assume to be spin-independent and the same for both 
leads. Since Ht conserves spin we have [H, S 2 ot ] — 0, and 
thus the two electrons from a given Cooper pair singlet 
which have tunneled to the lead(s) remain in the singlet 
state. 

Current of two electrons tunneling into different leads. - 
We use a T-matrix approach |2^| to calculate tunneling 
currents. At zero temperature the current I\ for tunnel- 
ing of two electrons coming from the same Cooper pair 
into different leads is given to lowest order in to by 



00 00 00 



h = 2e J dt j ' dt> j dt " e- ,,(t ' +t " 



)+i(2t-t' -t")fj. 



, -00 



x(H Tm (t ~ t")H Tm ,(t)H Tn (t')H Tn ,(0)), (3) 

where 77 — ^ + , and the expectation value is to be taken 
in the groundstate of the unperturbed system. The phys- 
ical interpretation of Eq. J3J) is a hopping process of two 
electrons with opposite spins from two spatial points ri 
and V2 of the SC to the two leads 1,2, thereby remov- 
ing a Cooper pair in the SC, and back again. The delay 
times between the two tunneling processes of the elec- 
trons within a pair is given by t' and t", resp., whereas 
the time between destroying and creating a Cooper pair 
is given by t. This process is contained in the correlation 
function 



£ ( H Tm<t ~ *") H Tm> (*) H T n(t') H Tn , (0)} 

n ? «' 

m 5* m' 

= M 4 £ {Gn«{t-t")G m ,- a {t-t')T nma {t')K ma {t'^ 

- G m ^ a (t -H- t") G na {t) T nma {t') K^-At") (e^C-*'"'") e -<*-C))< e ^»W e ~^)) . (4) 
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The lead Green's functions are G niJ (t) = (i/j na (t)ip^ a (0)} ~ (i/i/2)/it, with 2/j being the 
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DOS per volume at the Fermi level /i; of the 
leads. The anomalous Green's function of 

the SC is T nma {t) = (*_ CT (r m ,i)* ff (r„,0)) = 



(sgn(er)/Vs) J2k u kVk exp(-iE k t 



Sr) with 



Sr = ri — r 2 , and Vs is the volume of the 
SC. The bath correlator can be expressed as 
(exp(i(f> n (t))exp(-i<f> n (0))) = exp[J(i)] with 
J(t) = 2 f™(du/uj){ReZT(uj)/RQ)(exp(-iu}t) - 1). Here 
we introduced the total impedance Zp = (iu;C + R^ 1 )^ 1 , 
with a purely Ohmic lead impedance Z n (uj) = R, which 
we assume to be the same for both tunnel-junctions and 
leads. For small times, u> R \t\ <C 1, we can approximate 
J(t) ~ —iE c t where E c = e 2 /2C is the charging energy 
and uj r — l/RC is the bath frequency cut-off which 
is the inverse classical charge relaxation time r c / of an 
RC-circuit. For the long-time behavior, uj R \t\ 3> 1, we 

get J(t) (2/g)[ln(iu> R t) + 7] with 7 = 0.5772 the 

Euler number and g = Rq / R is the dimensionless lead 
conductance which determines the power-law decay of 
the bath correlator at long times. 

We first consider the low bias regime fi A,uj r . In 
this limit the delay times t' and t" ~ 1/A can be 

neglected compared to t ~ in all correlators in (@J 
and the bath correlators are dominated by the long-time 
behavior of J(t). We then obtain for the current 
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r(4/ 5 + 2) 



2/A 



4/9 



(5) 



The geometrical factor coming from the spa- 
tial correlation of a Cooper pair is Fd = z{8r) = 
[sm(kpSr)/kFSr] exp(— Sr/n^) with Sr = \Sr\. The 
exponential decay of the correlation sets in on the length 
scale of the coherence length £. It is on the order of mi- 
crometers for usual s-wave materials and can be assumed 
to be larger than Sr which could reach nanometers. 
More severe is the power-law decay oc l/(kpSr) 2 with 
kp the Fermi wavenumber in the SC. This power-law 
is sensitive to the effective dimensionality d [27j of the 
SC with weaker decay in lower dimensions. Indeed, in 
two dimensions and for kpSr 3> 1, but still Sr < £ we 
get F^_ 2 cx l/(kpSr) and in one dimension there is no 
power-law decay as a function of kpSr. In JSJ) we intro- 
duced the Gamma function T(x) and the dimensionless 
tunnel-conductance T = 7ri/gi/;|io| 2 with v$ being the 
DOS per volume of the SC at the Fermi level The 
result shows the well known power-law decay at low bias 
/i characteristic of dynamical CB [2^]. The exponent 
A/g in JSJ is two times the power for single electron 
tunneling via one junction. This is so, because the two 
tunneling events are not correlated since each electron 
tunnels to a different lead and the charge relaxation 
process for each circuit is independent. 
We consider now the large bias regime A,/i 3> lur. In 
the regime A, \fi — E c \ ^> ujr we can use the short time 
expansion for J(i) in As long as \fi — E c \ <C A we 
can again neglect the delay times t' and t" compared to t 
in all correlation functions in Q and obtain the current 



Ii in the large bias limit and up to small contributions 
~ evT 2 Fj{Sr)u R [0{u R /n) + 0(u> R /\n - E c \)} 



h = enr 2 F%(Sr)0(ii - E e )(ji - E c ). 



(6) 



This shows the development of a gap in I\ for fi < E c 
and R — ► 00 which is a striking feature of the dynamical 
CB. 

Current of two electrons tunneling into the same lead.- 
We turn now to the calculation of the current I2 car- 
ried by spin-entangled electrons that tunnel both into 
the same lead either 1 or 2. The current formula for I2 
is given by © but with n = n' and m = w! = n, and 
we assume that the two electrons tunnel off the SC from 
the same point and therefore Sr — here. Since both 
electrons tunnel into the same lead the bath correlation 
functions do not separate anymore as was the case in Q . 
Instead we have to look at the full 4-point correlator 



/ e #„(t-t") e #»(«) e -^«(*') e -#»(o)\ 
e J(t-t'~t")+j(t-t')+j(t-t")+j(t) 
~ e J(t')+J(-t") ' 



(7) 



The lead correlators again factorize into a product of two 
single-particle Green's functions since they are assumed 
to be Fermi liquids and in addition there appear no spin 
correlations due to tunneling of two electrons with oppo- 
site spins. 

We first consider the low bias regime \i ojr, A. Here 
again we can assume that \t\ is large compared to the 
delay times t' and t" , but it turns out to be crucial 
to distinguish carefully between A 3> ur and A 
lur. We first treat the case A 3> ojr and approximate 
exp[-(J(t') + J{-t"))] ~ exp[-iE c (t" - t')} in 0. In 
this limit and for A > E c the current I2 becomes 



9 (4A/tt) 2 „ 
e-KfjbT A „„ arctan 



El 



E r 



E r 
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r(8/ 3 + 2) \ojrJ 



8/9 



(8) 



The exponent 8/g in (JHJ we would also get in a first- 
order tunneling event if the operator e -11 ^" is replaced 
by e _l2< ^™ in (J2J) which changes the charge of the tunnel 
junction capacitor n by 2e. In addition to this double 
charging effect we find from (JHJ that an enhancement of 
E c gives not only rise to a suppression of I2 via the term 
(2/Vwr) 8/9 = {2^/gE c ) 8 /a but also to an increase due 
to the A-dependent prefactor. This enhancement can be 
interpreted as a relaxation of the charge imbalance cre- 
ated by the first electron tunneling event at small times, 
much smaller than the classical relaxation time t c ; . The 
result © is valid if >/(A - E C )/(A + E c ) > y/u) R /A. 
In the other limit where A <C u> R , e.g. for small R, we can 
assume that to R t' and oj R t" 3> 1 and therefore we approx- 
imate exp[-(J(f)+J(-t"))] ~ exp(-47/ 5 )w* /9 (i'i") 2/9 . 
In this limit we obtain 



(9) 
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with = (2e-T) 4 /sr 4 (l/ 3 +l/2)/7r 2 r(8/ 5 + 2). Here 
the relative suppression of the current I 2 compared to I\ 
is given essentially by (2/i/A) 4 ^ 9 and not by (2/j,/lor) 4 ^ 9 
as in the case of an infinite A. This is because the vir- 
tual state with a quasiparticle in the SC can last much 
longer than the classical relaxation time t c i, and, as a 
consequence, the power law suppression of the current is 
weakened since A <C lor here. To our knowledge, the 
result © was not discussed in the literature so far [30| . 
but similar results are obtained when SCs are coupled to 
Luttinger liquids Q . It is important to note that a large 
gap A is therefore crucial to suppress I 2 . 
In the large voltage regime A, [i 3> lor, we expect a 
Coulomb gap due to a charge q — 2e. Indeed, in the pa- 
rameter range \/i—2E c \ 3> lor and A 3> \fi—E c \ we obtain 
I 2 again up to small contributions ~ eirT 2 L0R{O(u>R/ /i) + 
0{co R /\fi-2E c \)} 

I 2 = eTTT 2 0{n-2E c ){n-2E c ). (10) 

This shows that I 2 is small (oc lo r /\^jl—2E c \) in the regime 
E c < (i < 2E C , whereas h is finite (oc F%(5r)(fi — E c )). 
Discussion and conclusions. - We now give numerical val- 
ues for the current magnitudes and efficiencies of our en- 
t angler. We first discuss the low bias regime (1 <C A,w_r. 
In Fig. 2 we show the ratio I2/I1 (efficiency of entan- 
gler) and I\ for A ^> E c ,ojr as a function of 4/g for 
realistic system parameters (see figure caption). The 
plots show that a very efficient entangler can be ex- 
pected for lead resistances on the order of R ~ Rq. 

The total current is then on the order of I\ ~ 10 fA. In 
the large bias regime \x » lor and for E c < fi < 2E C 
we obtain I 2 /h oc (k F Sr) d ~ 1 LO R /(2E c - fj,)(fi - E c ), 
where we assume that 2E C — fi and /1 — E c 3> lor. For 



fi ~ 1.5E C and using to r = gE c /-K we obtain approx- 
imately I 2 /Ii oc (k F 5r) d ~ 1 g 2 . To have I 2 /I\ < 1 we 




FIG. 2: Current ratio I2/I1 (entangler efficiency) and current 
Ii in the low bias regime, fi <C A, lor and A S> E c ,uor, as a 
function of 4/g = 4R/Rq. We have chosen realistic parame- 
ters: E c = 0.1 meV, k F Sr = 10, T = 0.1. The left plot is for 
/i = 5 /ieV and the right one for /1 = 15 /j,eV. In the case of a 
2d SC, Ji and I1/I2 can be multiplied by 10. 

have to achieve g 2 < 0.01 for d = 3, and g 2 < 0.1 
for d — 2. Such small values of g have been pro- 
duced approximately in Cr leads p"rj . For I\ we obtain 
h ~ e{k F 5r) 1 - d {^i - E C )T 2 ~ e{k F 5r) 1 - d E c T 2 ~ 2.5 pA 
for d = 3 and for the same parameters as used in Fig. 2. 
This shows that here I\ is much larger than for low bias 
voltages, but to have an efficient entangler very high lead 
resistances on the order R ~ 10i?Q should be used. Our 
discussion shows that it should be possible to implement 
the proposed device within state of the art techniques. 
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